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Abstract 

In a recent paper, we proposed a new class of supersymmetric SO{lQ) models for neutrino masses 
where the TeV scale electroweak symmetry is SU{2)l x SU{2)r x U{1)b-l making the associated 
gauge bosons Wr and Z' accessible at the Large Hadron Collider. We showed that there exists a 
domain of Yukawa coupling parameters and symmetry breaking patterns which give an excellent 
fit to all fermion masses including neutrinos. In this sequel, we discuss an alternative Yukawa 
pattern which also gives good fermion mass fit and then study the predictions of both models for 
proton lifetime. Consistency with current experimental lower limits on proton life time require the 
squark masses of first two generations to be larger than ~ 1.2 TeV. We also discuss how one can 
have simultaneous breaking of both SU{2)fi x U{1)b-l and standard electroweak symmetries via 
radiative corrections. 
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I. INTRODUCTION 



The nature of TeV scale new physics beyond the standard model (SM) is a question 
of enormous interest as the Large Hadron Collider (LHC) is poised to collect data in this 
energy range. Clearly, supersymmetry (especially the minimal supersymmetric extension of 
the standard model (MSSM)) is one of the prime candidates for this new physics since it not 
only solves the gauge hierarchy problem, but also has a number of attractive features such 
as the unification of gauge couplings at a high scale, a potential dark matter candidate, etc.. 
An interesting question along these lines has always been to see if any other new physics can 
co-exist with TeV scale supersymmetry without conflicting with coupling unification and 
dark matter, thereby broadening the scope of LHC physics search. 

A particularly appealing possibility is that weak interactions conserve parity asymptot- 
ically jll with the associated gauge group being SU{2)l x SU{2)ji x U{1)b-l so that the 
resulting gauge bosons Wr and Z' are at the TeV scale co-existing with supersymmetry. The 
case for SU{2)l x SU{2)r x U{1)b~l becomes more compelling when the SM or MSSM are 
extended to understand small neutrino masses via the seesaw mechanism {2!. As a generic 
possibility, this scenario is quite consistent with current low energy observations. Whether 
a TeV Scale SU{2)fi symmetry is compatible with supersymmetric couplirig unification has 
been extensively investigated in literature ^ With a few exceptions \4], it seems very 
hard to reconcile this possibility with the observed value of sin"^ 9w In a recent paper |5|, 
we pointed out a new supersymmetric 5*0(10) scenario where the presence of a vector like 
electroweak singlet and color triplet Higgs multiplet (which is part of the 45 representation 
in 5*0(10)) in addition to two bidoublets and two right handed doublets of the left-right 
electroweak group at the TeV scale leads to gauge coupling unification with TeV scale right 
handed Wr and Z' bosons . This model is different from other such scenarios considered in 



the literature j4| in that quark masses and mixingarise in a simple rnanner. The neutrino 
masses arise out of an inverse seesaw mechanism 6| and was shown [5] to have interesting 
phenomenological consequences like leptonic non-unitarity, leptonic CP-violation, lepton 
flavor violation, etc. which may be testable in near future. This fit to the fermion masses 
defines one class of 50(10) models with TeV scale which we call model (A). 

In this paper, several new results for these 5*0(10) models are presented: (i) we present 
an alternative fit to fermion masses, which we call model (B); (ii) we discuss the constraints 



of proton decay for both fermion mass fits - the one in Ref. [sl and the new one discussed 
in this paper; (iii) we also show how both B — L and electroweak symmetries can be broken 
radiatively in these models. 

Strength of proton decay has been studied extensively in the context of many supersym- 
metric grand unified theories (SUSY GUTs) (see Ref. {3] for recent reviews). Although there 
is no evidence for proton decay till now, current experimental lower bounds on the partial 
lifetimes of various proton decay modes tend to put severe constraints on these models e.g. 
they have now ruled out the simplest versions of SUSY SU (5) and suggest possible modifi- 

8| . They also constrain the choices of Higgs multiplets that can be 



cations of such models 



used for model building with 5*0(10) group jo]. 

In the models we are discussing here, due to the fact that all the Yukawa couplings 
responsible for proton decay are constrained by the fermion mass fits, it is possible to 
estimate the partial life times for the various modes as functions of the squark masses and 
for reasonable squark masses of the first two generations, and for model (A), we get upper 
bounds on various proton decay channels. There are no such bounds in the second case 
(model (B)). We find that within a reasonable set of assumptions, all our predicted upper 
bounds for model (A) are consistent with the current experimental bounds and some of the 
modes may be accessible to the next generation proton decay experiments with megaton 
size detectors. 

We also discuss the constraints imposed by radiative breaking of both SU{2)ji x U{\)b-l 
and the SM gauge symmetries via radiative corrections. The idea is to start with soft mass 
squares at the Planck or GUT scale and extrapolate the masses to the weak scale to see if 
the SU{2)fi X U{1)b-l symmetry breaks at the TeV scale. We then note that this breaking 
introduces via D-terms a breaking of the SM gauge symmetry to f/(l)em- 

We also discuss the generalization of this model to include i?-parity breaking and its 
implications on proton decay. 

This paper is organized as follows: in Sec. II, we review the basic structure of our model 
and the gauge symmetry breaking. In Sec. Ill, we review the fermion mass fit for model (A) 
already discussed in Ref. jsl. In Sec. IV, we present a new fermion mass fit and define it as 
model (B). Sec. V describes the radiative electroweak symmetry breaking (EWSB) in this 
type of models. In Sec. VI, we discuss the proton decay in both these models. In Sec. VII, 
we comment on the effect of i2-parity breaking terms in the superpotential on proton decay. 



The results are summarized in Sec. VIII. In Appendix A, we present the renormalization 
group equations (RGEs) for soft SUSY-breaking masses in our supersymmetric left-right 
(SUSYLR) model. In Appendix B, we derive the anomalous dimensions of the dimension-5 
proton decay operators in our model. In Appendix C, we list the hadronic form factors used 
in our proton decay calculations. 



II. A BRIEF OVERVIEW OF THE MODEL 



As in the usual 5*0(10) models, the three generations of quark and lepton fields are 
assigned to three 16 dim. spinor representations. In addition, we add three 5*0(10) singlet 
matter fields to implement the inverse seesaw mechanism. The B — L gauge symmetry is 
broken at the TeV scale by 16-Higgs fields (denoted by iPh)i whereas the rest of the gauge 
symmetry is broken at ~ 10^^ GeV by 54 and 45- fields (denoted by E and Aa respectively). 
We require two 45-IIiggs fields (a = 1,2), one for symmetry breaking and the other to give 
rise to the vector-like color triplets at the TeV scale. The SM symmetry is broken by two 
10-Higgs fields (denoted by Hg)- We note that the field content of our model is found in 
many string models after compactification e.g. fermionic compactification models jl^ and 
it may therefore be easier to embed this GUT model into strings. 

The distinguishing feature of our model is that the GUT symmetry breaks down to the 
left-right symmetric gauge group 5f/(3)c x SU{2)l x SU{2)ji x U{1)b-l without parity (Im- 
parity). The D-parity is broken at the GUT scale by the vacuum expectation value (VEV) of 
the 45-Higgs field. A consequence of /^-parity breaking is that only the right-handed (RH) 
doublets from 16-Higgs fields survive below the GUT scale. An interesting feature of this 



class of models {s] is that if we have two RH Higgs fields [x'^,x'^ (1, 1, 2, ±1)] , two bi-doublet 
fields [$(1,2,2,0)] (all color singlets) and a vector-like color triplet but SU(2)l x SU{2)r 
singlet field [6 (^3, 1, 1, |^+c.c.] at the TeV scale, the gauge couplings unify around 10^^ GeV. 
The bidoublet fields arise from 10-Higgs at the GUT scale and the vector-like color triplet 
fields arise from the 45-Higgs field. This is therefore a new class of 50(10) SUSY-GUT 
theories with TeV scale Wr and Z' bosons which can be accessible at the LHC. 
We consider the symmetry breaking chain 

^O(IO) ^ 3,2l2r1b-l ^ 3,2^1y(MSSM) 3,2Uy(SM) ^ 3,1q (1) 

where, as an example of our notation, 3^ means SU (3)c. As shown in Appendix A of Ref. 



for consistency, we need at least two 45 and one 54 representations of the Higgs fields to 
break the SO{10) gauge group into the SUSYLR gauge group, SU{3)c x SU{2)l x SU{2)r x 
U{1)b~l, at the scale Mq ~ 4 x 10^® GeV. Note that to have realistic fermion masses and 
mixing, we need at least two SU (2) bi-doublets of the 10 Higgs representation to break the 
SU{2)l X U{1)y gauge group of the SM to U{1)q at the weak scale Mz- With this minimal 
set of Higgs fields, we were able to attain not only gauge coupling unification but also the 
desired fermion masses and mixing at the GUT scale . Incidentally, since our gauge group 
above TeV scale is different from MSSM, we needed to extrapolate fermion masses using the 
left-right group (see Appendix B of Ref. 5|]) which has certain distinguishing features in the 
running behavior, in contrast to the MSSM gauge group. 
The superpotential for the model consists of several parts: 

W = WsB+ +W' (2) 

where Wsb is responsible for 50(10) GUT symmetry breaking, doublet triplet splitting and 
the remnant sub-GUT scale multiplets; is the Yukawa superpotential responsible for 
fermion masses and mixing; W involves the i?-parity violating terms. When we impose an 
additional matter parity symmetry under which ip^ — ^ —i^a, — )• — 5'q,, and all other fields 
even, as was assumed in Ref. 5(], we get W = 0, i.e. all -R-parit v v iolating terms are absent 
in the superpotential and the model has a stable dark matter ll|. We discuss the effects 
of nonzero W in a subsequent section where we show that even after including arbitrary 
i?-parity violating terms (i.e. giving up matter parity assumption), the model does satisfy 
proton life time bounds since W conserves baryon number and after B — L breaking leads 
to a highly suppressed amplitude for proton decay. This feature is characteristic only of 
5*0(10) models with low B — L breaking. 
The Yukawa superpotential is given by 

= Ki.lQ^lQ.lQu^ + ^16a6,10^^„45H45'^ (3) 

where the first term is the usual Yukawa coupling term, while the second term is a higher- 
dimensional term whose completely antisymmetric combination acts as an effective 126if 
operator, thus giving rise to a realistic fermion mass spectrum at the GUT scale. We define 
this as our model (A). 

The superpotential Wsb was discussed in detail in Ref. jsl where it was noted that the 
following components of the 54 and 45 Higgs fields acquire VEV and leave the left-right 



subgroup unbroken: 



A- f 3 3 3 3 

(54) = diag la, a, a, a, a, a, "2*^' "2*^' ^2*^' "2*^ 

(45)12 = (45)34 = (45)56 = h. (4) 



III. FERMION MASSES IN MODEL (A) 

The model discussed in Ref. 5| is defined by the VEV pattern of tlie bi-doublets: 



(<f 1) = , {^2) = (5) 

\ y \Q K^j 

We define the ratio of the VEVs as tan/3 = — as in MSSM. Then the fermion mass matrices 
at the GUT-scale are given by 

Mu = ~K + /, 
Md = hd + /, 
M, = hd- 3/, 
M,^ = K- 3/. (6) 

where in the notation of Ref. [sj, hu4 = i^u,dhu,d- The contribution from the effective 
operator is assumed to be the same for both up and down sectors, i.e. / = Hufu = f^dfd'-, as 
a resuh, we have the relation fd = fu tan /3. Also note that the factor —3 between the quark 
and lepton sector is due to the same 126 operator. Using the renormalization group analysis 
for the fermion masses and mixing in the SUSYLR model (see Appendix B of Ref. [5]), we 
obtain the GUT-scale fermion masses starting from the experimentally known values at the 
weak scale. Using these mass values, we obtain a fit for the coupling matrices at the GUT 
scale defined in Eq. ([6]). Here we give the results in a down quark mass diagonal basis for 
two cases: 

(a) tan/^MSSM = 10: In this case, the GUT-scale values of the charged fermion masses are 
found to be 

ruu = 0.0017 GeV, m, = 0.1908 GeV, rrit = 77.7 GeV, 
nid = 0.0013 GeV, m, = 0.0263 GeV, ruf, = 1.7001 GeV, 
me = 0.0004 GeV, = 0.0910 GeV, = 1.7061 GeV (7) 



and tan /3gut = 7. Note that the GUT-scale fermion masses quoted here are shghtly different 
from those given in Ref. jsl because, in this case, we have set the 5$$ couphng = (of 
Ref. jsl) assuming i?-parity conservation. With these mass eigenvalues, we find a fit for the 
GUT-scale couplings of the form: 

U = diag (l.26 X 10-^ -0.0001, -9.48 x 10^^) , fd = /.tan/^cuT, 
hd = diag (4.86 x 10"^ 0.0019, 0.0752) , 

^ 7.46 X 10-5 0.0002 - 6.51 x lO'^i 0.0002 - 0.0028z ^ 



0.0002 + 6.51 X 10- 



0.0015 



0.0118 + 1.26 X 10- 



(8) 



y 0.0002 + 0.0028Z 0.0118 - 1.26 X lO-^i 0.4908 

Note that for simplicity we have chosen the /-couplings to be diagonal. Our fit does not 
allow the off-diagonal components to be too different from zero. 



(b) tan/^MssM = 30: In this case, the GUT-scale values of the charged fermion masses 
are found to be 

rriu = 0.0121 GeV, = 0.3269 GeV, = 120.53 GeV, 
rud = 0.0014 GeV, = 0.0277 GeV, rrib = 2.7958 GeV, 



me = 0.0006 GeV, = 0.1266 GeV, = 2.7737 GeV 



(9) 



hu 



and tan/^GUT = 20. With these mass eigenvalues, we obtain a fit for the couplings of the 
following form: 

fu = diag (l.5 X 10-^ -0.0002,4.2 x 10-^) , fd = /^tan/^cuT, 
hd = diag (0.0002,0.0078,0.4163), 

0.0002 0.0003 - O.OOOli -0.0008 - 0.0081i ^ 

0.0002 + 0.0001Z 0.0029 0.0144 + 0.0002i (10) 

-0.0008 + 0.008n 0.0144 - 0.0002i 0.9145 j 

We note that in this model, larger values of tan/3 (> 30) are not allowed. This can be seen 

n 

analytically from the form of the RGEs given in Appendix B of Ref. [5| where it is clear 
that the up-quark sector masses will increase rapidly at high energies for large tan /3 and the 
same effect is induced in the down-quark sector which makes the Yukawa terms dominant 
over the gauge terms. This makes all the quark masses to run up to unacceptably large 



V 



values at the GUT-scale. We believe this is a general feature of low-scale SUSYLR models, 
in contrast to MSSM case 12 1 . 



IV. A NEW FERMION MASS FIT: MODEL (B) 

In this section, we consider an alternative mass fit within the 5*0(10) models with low 
scale B — L. It follows from a recent ansatz [l^ that in generic 5*0(10) models which do not 



use type 
the idea 



seesaw to fit neutrino masses, an alternative fit to fermion masses is possible using 



3 has a rank one 10-Higgs Yukawa coupling matrix which dominates 

the fermion masses while other couplings introduce small corrections; the third generation 
masses arise from the dominant rank one coupling matrix with smaller 126 and second 10 
couplings generating the CKM mixing as well as the second and the first generation fermion 
masses. This idea can be applied to our case since, the neutrino mass is given by the inverse 
seesaw formula which involves an additional matrix /i. The main difference of model (B) 
as compared to model (A) resides in the VEV pattern of the two Higgs bidoublets i.e. in 
model (B), we have 



^ f^d 

^ 



< 
k' 



(11) 



with f„k/-\/2 = J k'I + + + i^'d- Also we must have — 7^ ^ in order to get right 
fermion mixing pattern. In the limit k„ ^ k^, the RG analysis of model (A) can be applied 
to this case to generate fermion masses at the GUT scale as well as the symmetry breaking 
pattern via radiative corrections. 

The resulting fermion mass formulae in terms of the appropriately redefined Yukawa 
couplings are given as follows jl^ : 

Md = n{h + f + h'), 
Mi = n{h-^f + c,h'), 
M,^ = h - 3f + c,h' (12) 

where 



Hd Hd ' 



Ti = — , r2 = ra = -. 



(13) 



As in the case of model (A), the / couphng above represents the effective 126 couphng 
arising from the 1IJ1IJA1A2H2 term in the supcrpotential and h' arises from a coupling of the 
form 1IJ1IJH2X (with a nonzero VEV for the additional singlet field X). Note that if there 
is an additional Z2 symmetry under which H2,A2,X are odd and all other fields are even, 
one can have a supcrpotential with only the h, f, h' type contributions as given above, to 
the fermion mass formulae. In our case with two Higgs bi-doublets, Cg — 1 and Ci, — r^. 
With the GUT-scale mass eigenvalues obtained earher, we obtain a fit for these couphngs 
as follows: 
(a) tan/^MSSM = 10: 

Ku = 173.2 GeV, n = 0.0218, ra = 0.14, 
/i = diag (0,0,0.45), 

^ -0.0006 0.0019 ^ 



/ 



-0.0006 0.0115 0.0101 
0.0019 0.0101 0.0001 



h' 







-0.0022 0.0005 ^ 
0.0022 0.0181 
-0.0005 -0.0181 y 



(b) tan/^MSSM = 30: 



(14) 



Ku = 172.4 GeV, n = 0.0231, rs = 0.21, 
/i = diag (0,0,0.70), 

^ -0.0016 0.0062^ 
/= -0.0016 0.0140 0.0111 (15) 

^ 0.0062 0.0111 0.0019 y 

and h' same as in case (a). It may be noted here that in both the cases, all the fermion mass 
values predicted using the couplings above agree with those obtained from the RGEs within 
the experimental uncertainty, the only exception being the up-quark mass in case (a), where 
the our predicted value is about 4 times larger. Note however that in our discussion, we 
have not included contributions from threshold corrections or higher dimensional operators. 
Those contributions can generally be of order MeVs when their couplings are chosen ap- 
propriately, in which case, they will not affect the second and third generation masses but 
could easily bring the up quark mass into agreement with RGE predictions. 



With the Yukawa couphngs completely fixed in our model, we can analyze the predictions 
for the proton decay rate. But before doing so, we discuss the details of the electroweak sym- 



metry breaking in this model which was not done in the original paper js]. This discussion 
applies to both models (A) and (B). 

V. SYMMETRY BREAKING BY RADIATIVE CORRECTIONS 

In this section, we propose a way to break both the SU{2)ff x U{1)b-l as well as the SM 
symmetry via radiative corrections from renormalization group extrapolation of the scalar 
Higgs masses from the GUT to TeV scale. As is well known, the large top quark coupling 



enables us to achieve a similar goal i.e. radiative EWSB in the case of MSSM [15|. The 
simple generalization of that procedure cannot work in our model since the bidoublet Higgs 
of LR models contains both the H^^d components of MSSM, and as a result, large top quark 
coupling will necessarily turn both their masses negative and this is known not to give a 
stable vacuum. 

Our proposal is that we use a domain of parameter space for the soft SUSY-breaking 
mass squares for the RH Higgs doublets x'^ where the mass square of one of them 

turns negative, by RG running to the TeV scale due to the L'^x'^S Yukawa coupling being 
large. This leads to a breaking of the SU{2) and B — L symmetry. The mass square of 
the x'^ remains positive throughout but it acquires an induced VEV. The differences in their 
VEVs, via the D-term, can make the mass square of the Hu field negative while keeping the 
mass square of positive as in the case of MSSM, thereby also giving rise to the EWSB. 
The main point is that both symmetry breakings owe their origin to one radiative correction. 

In order to show that it is indeed possible to achieve negative mass square for one of the 
RH Higgs doublets while keeping all other soft mass squares positive, we need to examine 
the RG running of all the soft mass parameters from the GUT to TeV scale. In this regime, 



the mode 



given by [16 1 



is SUSYLR for which the superpotential and soft SUSY-breaking Lagrangian are 



+iM^.r + /i$„,^"Tr [^r.^.r^) + M$„,Tr [^ir^^.r^ 
Ar^ap^iS^gPgi + ^MfS'^S^, (16) 
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Aoft = -\ (Mfg~g + M^lWlWl + M^rWrWr + M^BB + h.c. 



'T 



V 



1 



(17) 



where we have suppressed the generational and SU (2) indices, and a, 6 = 1, 2 (for two bidou- 
blets), p,q = 1,2 (for two SU{2)r doublets) and q;,/3,7 = 1,2,3 (for three gauge singlets). 
Note that we do not have any x-term in these expressions as there is no SU{2)l Higgs 
doublet in our model. Also we have an additional term in the superpotential (the SL^x'^ 



term) and a corresponding trilinear term in t 
as compared to the expressions given in Ref. 



le soft breaking Lagrangian (the SL'^x'^ term) 



16| ; this additional term in the superpotential 



is required for the inverse seesaw mechanism to work. Moreover, if we assume i?-parity 
conservation, then the Sx'^x'^ and 5*$$ terms are not allowed in the superpotential and also 
in the soft-breaking Lagrangian, i.e. the couplings /i^c and fiq, as well as Yabc in Eq. (ITBI) and 
the corresponding terms in Eq. (fTTI) are set to zero and /xlc is the only non-zero coupling in 
Eq. f|T6l) which can be fixed by requiring b — r unification at the GUT-scale. In this section, 
we work with this assumption; the effects of R- parity breaking will be discussed later. 

Now we analyze the RG evolution of the gaugino and soft mass parameters from GUT 
to TeV scale. It is well known that in SUSY GUTs, the /3-function for the gaugino mass is 
proportional to the /3-function for the corresponding gauge coupling. Explicitly, the RGEs 
for the gaugino mass parameters are given by 



dMi 



2bi 
IGvr^ 



(18) 



where the /3-function coefficients in our SUSYLR model are [5| hi = (13,2,4,-2), corre- 
sponding to i = 1b-l, "^l, 2ij, 3c respectively. This implies that the three gaugino masses, 
like the three gauge couplings, must unify at fi = Mgut- In order to solve Eq. ( fT8l) . we 
adopt the universality hypothesis at the GUT scale (as in typical mSUGRA type models) 



Ml = M2L = M2R = Ms = mi/2. 



(19) 
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together with the initial condition 



9I = alh = 9m = 9l = 4vraGUT, (20) 
where Mgut ^ 4x 10^^ GeV and 

'^GUT — 20.3 in our model {sl- Using these initial conditions, 
we can obtain the running masses for the gauginos at TeV scale, starting with a given value 
mi/2 at the GUT scale, as shown in Fig. [1] for a typical value of mi/2 = 200 GeV. The value 
of M3 increases, since it has a negative /3-function, while the other gaugino masses decrease 
as we go down the energy scale. Thus the gluino is much heavier than other gauginos at the 
weak scale. 
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FIG. 1. RG evolution of gaugino masses from GUT to TeV scale for mi/2 = 200 GeV. 

The one-loop RGEs for the soft SUSY-breaking mass parameters are given in Appendix 
A. As initial conditions, we assume universality and reality of the soft fermion and Higgs 
masses at the GUT-scale, i.e. 

^X" = = "^0' = ^l^ab, (21) 

whereas a different scale is assumed for the soft singlet scalar mass: 

{ml)^^=m'^ Va,/3 = 1,2,3. (22) 

In principle, we can choose a different mass scale for the Higgs bidoublets and even different 
generations of fermions as well. The only constraint due to the 5*0(10) symmetry requires 

12 



us to have the same mass for each generation of fermions. Note that all the off-diagonal soft 
SUSY breaking scalar masses have been set to zero. The inter-generation mixing at the low 
energy scale then occurs only via the superpotential Yukawa couplings. With these initial 
conditions, we solve the coupled RGEs for the soft masses given in Appendix A, along with 
the Yukawa RGEs given in Ref. to get the running soft masses at the low scale. We find 
that it is indeed possible to find a parameter space such that m|c < (for S'[/(2)r breaking) 
and < (for EWSB) while keeping all other mass squares positive. Fig. [2] illustrates 
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FIG. 2. Evolution of the scalar mass parameters for mi/2 = 200 GeV, niQ = 1.20 TeV and 
uiq = 1.27 TeV. For the scalar masses, we actually plot sign(m^) • y^|rn?|, so that the negative 
values on the curves correspond to negative values of m?. 

such a scenario for the choice mi/2 = 200 GeV, itlq = 1.2 TeV and rriQ = 1.27 TeV. We 
have chosen the SL'^x'^ coupling = 0.7 to achieve a realistic fermion mass spectrum, 
and in particular, the h — t unification at the GUT scale. Note that the RH slepton masses 
evolve much more rapidly than their LH counterparts due to this large coupling /i^c The 
value of ttlq is chosen such that all the other eigenvalues (especially m|c and m|) remain 
positive at the TeV scale. Note that the low energy values of m|c and m| are of order (10 
GeV)^. However the physical masses of these particles also receive a contribution from the 
{yf) which pushes the masses upto a TeV scale. As far as the squark masses are concerned, 
they evolve more than the slepton masses due to the strong interaction loop contributions 
to their RGEs. The small intra-generational mass splitting is due to the differences in 
their electroweak interaction. We can see clearly that at the weak scale, the values of 
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m|c and are negative, thus triggering the SU{2)ii and electroweak symmetry breaking 
respectively. Note that we need not have both the bidoublet mass squares to be negative, as 
one negative value will induce the symmetry breaking via the cross terms of the type $i$2 
in the Lagrangian. 

We also verify that the low-energy values of the sfermion mass square matrices satisfy all 
the FCNC constraints , due to the smallness of the off-diagonal entries. As an example, 
we give the values here for the parameter values shown in Fig. 2: 



ml 



mlc 



( 



1.63 X 10^ 



-1.45 X 10^ - 8.64 X lO^i 
-4.79 X 10^ - 3.57 x lO^i 



1.58 X 10^ 
-1.45 X 10^ - 8.64 X lO^i 
-4.79 X 10^ - 3.57 x lO^i 

1.39 X 10^ 
-7.28 - 8.39 X lO^i 
-2.59 X 10^ - 3.45 x lO^i 

3.81 X 10^ 
-7.18 - 8.24 X lO^i 
-2.57 X 10^ - 3.42 x lO^i 



1.45 X 10^-^8.64 X lO^i - 


-4.79 X 10^ + 3.57 x lO^i 




1.63 X 10^ 


-2.31 X 10"^ + 1.68i 




-2.31 X 10^ - 1.68f 


6.51 X 10^ 
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1.45 X 10^-^8.64 X lO^i - 


-4.79 X 10^ 3.57 x lO^i 
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1.58 X 10^ 


-2.31 X 10^ 1.68i 




-2.31 X 10^ - 1.68i 


5.99 X 10^ 


/ 


-7.28 8.39 X lO^i 


-2.59 X 10^-^3.45 x 10 




1.39 X 10^ 


-1.25 X 10^-^7.45 X 10" 


-\ 


1.25 X 10*^ - 7.45 X 


8.66 X 10^ 




-7.18 -h 8.24 X lO^i 


-2.57 X 10^-^3.41 X 10 




3.81 X 10^ 


-1.24 X 10^-^7.75 X IQ- 




1.24 X 10^ - 7.75 X 


5.00 X 10^ 





GeV^ 



GeV^ 



GeV^. 



VI. PROTON DECAY 

In this section, we discuss the partial lifetimes of various proton decay channels. 

A. Proton decay operators 

In generic SUSY-GUTs, there exist three sources for proton decay: 
• D-type (dimension-6) operators that arise from exchange of gauge bosons: 

-\ f d^e d^e $t$$t$, (23) 
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which may be generated both by heavy gauge boson exchange and by heavy chiral 
(Higgs) superfield exchange. For a unification scale ^ 10^^ GeV, these contributions 
to proton decay are sufficiently small and well beyond the range of current experiments. 

F-type (dimension-5) operators that arise from the exchange of color triplet Higgsino 
fields in 10-Higgs fields as shown in Fig. [3t^a): 

1 



M, 



G 



d'd $$$$ 



(24) 



where $'s are used to denote quark and lepton doublets. In the component language, 
they give rise to dimension-5 operators of the form {QQ){QL) and {QL){QQ). As 
these operators involve squark and slepton fields, they cannot induce proton decay in 
the lowest-order. Proton decay occurs by converting the squark and slepton legs into 
quarks and leptons by exchanging a gaugino, as shown in the box diagram of Fig. 
Hb). 

Another class of dimension-5 operators arising from i?-parity breaking Planck sup- 
pressed operators, which are absent when we assume i?-parity. We discuss them in 
Sec. VI and show that their effects are very small due to low B — L breaking scale. 
These are absent in models where 126 Higgs fields break B — L, but are present in 
our model. 

Q Qk 




Wh Wh 

X " 



(a) 




h: 



9: W, B 



(b) 



FIG. 3. (a) Supergraph giving rise to effective dimension-5 proton decay operators, and (b) Box 
diagram involving gaugino exchange that converts the dimension-5 operator of Fig. 3(a) into an 
effective four- Fermi operator that induces proton decay. 



There are two effective dimension-5 operators of LLLL type that involve only left-handed 
quark and lepton fields, given by Eq. (12^ and a corresponding RRRR type, both invariant 
under MSSM |l^. In super-space notation, these are explicitly given by 



(25) 
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On = l d'e 6"^^ iQ^)^^ (Q^)^^. m,, {L^)i 



(26) 



where a,/3,7 = 1,2,3 are SU{3)c color indices; a,b,c,d =1,2 are SU{2)i isospin indices; 
and i,j, k, I = 1,2,3 are generation indices. It is clear from the form of these operators that 
they break baryon number by one unit, but preserve the B — L symmetry, leading to the 
proton decay to a pseudoscalar and an anti-lepton. As argued in Ref. 13] for kinematical 



reasons and explicitly shown in Ref. 



20[ | for small to moderate tan /3 region of the SUSY 



parameter space, the RRRR contributions are at least an order of magnitude smaller than 
the LLLL contributions. We also verify this in our model, as shown later; for the time being 
therefore, we concentrate only on the LLLL operator. 

In component form, the effective superpotential due to the LLLL operator is explicitly 



given by 



2l| 



A_B=1 



[{Cijkl — Ckjil) UaidfijU^k^l — {Cijkl — Cikjl) Uaidfijd^kl^i] (27) 



where is the effective mass of the color triplet Higgs field belonging to the lOi;;- repre- 
sentation, and in our model, is of the order of the unification scale Mq (see Appendix A 
of Ref. 5|). This superpotential leads to the effective dimension-5 operators involving two 
fermions and two sfermions as shown in Fig. El^b), which lead to proton decay by four- Fermi 
interactions when "dressed" via the exchange of gauginos, namely gluinos, binos and winos. 
A typical diagram for the effective four-Fermi interaction induced by this dressing is shown 
in Fig. m 




FIG. 4. The effective four-Fermi interaction diagram induced by the gaugino dressing of the 
effective dimension-5 operator given by Fig. [3)[^b). 



The coefficients Cijki associated with the superpotential given by Eq. (127|) can be ex- 
pressed in terms of the products of the GUT-scale Yukawa couplings. For model (A), this 
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ij JUkl I J Uij ' "Uf^l 

h 



is given by 

Cijki = K.jKki + xihd^Md^^ + X2hu,^hd^i + X3hd^Mu^^ + - / 

+xi [hd^Jdki + fd,,hdi,i) + X2 [fu,,hd^i + K^Jdki) + X3 [hd^Juki + fd,jK 

+ ^ [fuijfu,,l + Xifdijdki + X2fuijfdki + Xsfd^jfuki) (28) 

while for model (B) this becomes 

Cijki = hijhki + xih'ijh'^i + X2hijh'f^i + x^h'^jkhi + ^ ^1 {Kjfki + fijKi) + X2hijfki + x^fijhki 
+\xifijfki (29) 

where the ratios of the 10^ color triplet Higgs masses and mixings and the factor 

I is the C-G coefficient for the 10 ■ 10 ■ 126 coupling. Note that there are only three 
mixing parameters as there are only four color triplet Higgses in the MSSM gauge group, 
corresponding to the two 10// fields in our model. As we are interested only in the upper 
bound for the partial lifetimes of various proton decay channels, we do not need to know 
the detailed form for the Xi parameters in terms of these masses and mixings. We just vary 
these parameters numerically to get the maximum value for the partial lifetimes. 

It can be shown that 22] in the limit of all squark masses being degenerate as in typical 
mSUGRA type models, the gluino and bino contributions to the dressing of the dimension-5 
operators vanish. This basically follows from the use of Fierz identity for the chiral two com- 
ponent spinors representing quarks and leptons. In realistic models, the FCNC constraints 
allow only very small deviations from universality of squark masses. Hence, these gluino and 
bino contributions are expected to be small compared to the wino contributions, and can 



be ignored altogether. The charged wino dressing diagrams have been evaluated earlier 



23|, 



and in the limit of degenerate squark masses, this leads to the effective Lagrangian 21] 

Ub=i = 2Ie''^\Ck,a - Q,ki)[ulkCdp,(f^fiui + ulfd.,kulfieil (30) 
where C denotes the charge-conjugation matrix and I is given by 



in Mj ' 



(31) 



m^y being the wino mass and Mj the sfermion mass. Using this expression and adding a 
similar contribution from the neutral wino exchange diagram, we can write down the total 
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contribution to various proton decay channels. This is summarized in Table [B We note 
that the proton decay operators with s-quark lead to i^'-meson final states whereas the ones 
without s lead to vr final states. As shown in Table HI the amplitude for non-strange quark 
final states will be Cabibbo-suppressed compared to the strange quark final states. It is also 
important to mention here that the total amplitude for final states involving neutrinos is 
the incoherent sum of the rates for all three neutrino states. This leads to large decay rates 
for p — 7- K~^T7 and p — t- 7r+z7 channels compared to the other decay channels due to the large 
Yukawa couplings of the third generation. 



Decay channel 


C-coefficient 


p — )■ K^Vi 


{Cll2l — Ci2ll) 


p K°e+ 


(C1121 — C1211) 


p K^fi+ 


(C1122 — C'1212) 


p — 7r"*"I7; 


sin6'c(C2iii - C1121) 


p — )• vr'^e"'" 


sin0c(C2iii — C1121) 


p —7- TT^ n'^ 


sin6'c(C2ii2 - C1122) 



TABLE I. The coefficients for various /S.B = 1 dimension-5 operators obtained from the effective 
Lagrangian to leading order. Here 6c is the Cabibbo angle (with sin^c" ~ 0.22) and the Cijki's are 
products of the Yukawa couplings, as defined in Eqs. (|28|) and ([29]) . 



Before proceeding to calculate the rate of proton decay induced by these LLLL type 
operators, let us estimate the contribution from the RRRR type operators in our model. 
The gluino dressing graphs do not contribute in the limit of universal sfermion masses by the 
same Fierz arguments as for the LLLL case. Moreover, since all superfields in the RRRR 
operator are SU(2)l singlets, there is no wino contribution to the leading order. Also the 
bino dressing generates an effective four-Fermi operator of the type e°''^''e'^^e^''ufjCd'^i^u'^^Cef 
which, in fiavor basis, is antisymmetric in the fiavor indices i and j, and hence in the 
mass basis, must involve a charm quark. Thus to leading order, the bino contribution 
also vanishes due to phase space constraints. Thus the only dominant contribution comes 
from the Higgsino exchange and the largest amplitude in this case, which comes from stop 



intermediate states, is estimated to be 



21| (using the Cijki values calculated later in our 
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model) 

gi323 ,, , - 4.0 X 10-^° (32) 

IGTr^'f 4k sm [5 cos p 

for tan/3 = 30, as compared to the LLLL contribution which is typically of order 



Cii23^ ~ 4.5 X 10-9 (33) 

As the RRRR contribution is proportional to ^.^^'''^^^^ which is ~ tan/3 for large /3, for 
smaller tan/3, this contribution is further suppressed. This justifies why we can ignore the 
RRRR contributions in the following calculation of proton decay rate. 



B. Proton decay rate 

In order to calculate the proton decay rate, we must extrapolate these dimension-5 op- 
erators defined at the GUT scale to the scale of ntp = 1 GeV. In our model, we can divide 
this whole energy range into three parts, following the breaking chain given by Eq. ([T]): 

(a) from the GUT scale Mq to the B — L breaking scale M/j (SUSYLR), 

(b) from Mr to the SUSY-breaking scale Ms (MSSM), and 

(c) from Ms to 1 GeV (SM). 



The values of these extrapolation factors are given in the literature [19|, l24j-|26| for both 
SM and MSSM, but not for the SUSYLR model. In this section, we derive these factors 
using the anomalous dimensions for the dimension-5 operators in our model calculated in 
Appendix B. We denote the overall extrapolation factor by Af,. We noted some discrepancies 
in the values of the anomalous dimensions quoted in different papers, but found that our 
results for the SM and MSSM cases agree with those given in Refs. Q, Q and quoted in 



Appendix E of Ref. 

We also need to include the QCD effects in going from three quarks to proton. As the 
low-energy hadrons are involved in the decay, this is a highly non- perturbative process, 
and it is difficult to calculate the exact form of the hadronic mixing matrix element for 
the process. Even though various QCD models have been constructed for the purpose, 
the estimates vary by a factor of (9(10) between the smallest and the largest As the 
partial width of the decay is proportional to the matrix element squared, the variation in the 
estimate of proton lifetime in different models will be (9(100). A different approach using 
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lattice QCD techniques gives more consistent results [28[. We use these recent results to 
estimate the chiral symmetry breaking effects which can be parametrized by two hadronic 
parameters D and F. Then the hadronic mixing matrix for the proton decay can be written 
as -j^f{F, D) where = (130.4 ± 0.04 ± 0.2) MeV 3 is the pion decay constant and 

= 0.0120(26) GeV^ 28| is a low-energy parameter of the SU (3) / baryon chiral Lagrangian 
with the baryon number violating interaction. The factors /(-F, D) for different final states 
are listed in Appendix C. 

Finally, combining all the factors discussed above, the proton decay rate for a given decay 
mode p — )■ Ml (M denotes the meson and / the lepton) is given by [21 1 

2 



Tp{Ml) 



1.6 X 10"^^ GeV 
X \C\^\A,\^\f{F,D 




(34) 



where the coefficients C are given in Table [H the hadronic factors f{F,D) are listed in 
Appendix C, and the extrapolation factors A^. are derived below. 



C. The extrapolation factors for the dimension-5 operator 



As noted in the previous section, we need to extrapolate the dimension-5 operators defined 
at the GUT scale to the scale of 1 GeV. In our model, this whole energy range is divided into 
three parts, with different running behavior for the gauge couplings. First, we have the SM 
sector from 1 GeV to the SUSY-breaking scale Ms in which we have the usual non-SUSY 

n 

enhancement factor [24] for the LLLL operator: 



A 



NS 



03(1 GeV) 



2/(ll-|"/) 



(35) 



a^iMs) 

where nj is the number of quark flavors below the energy scale of interest. Here we have 
neglected the effects of SU{2)l and f/(l)y couplings as they are much smaller compared to 
that of S'f/(3)c. In our model, as Ms = 300 GeV > rrit, the enhancement factor explicitly 
becomes 

-| 2/7 

= 1.49 (36) 



A 



NS 



"a3(l GeV)' 


2/9 




6/25 


a3(mb) 


6/23 


03 (mt) 






_ot^{mb)_ 
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using the values of as{fi) at /i = 1 GeV, rric and obtained by interpolating the renor 



30| and ai = rrit by the SM 



malization group equation for the effective QCD coupling 
running from /i = m^. 

Now above Ms, we have the usual MSSM till the B — L breaking scale Mr and then the 
SUSYLR model till the GUT scale Mq- The extrapolation factor in this case is given by 



MSSM /I SUSYLR 



/lb _ ^MbSM A, 



(37) 



where the corresponding factors in the two sectors are given by 



A 



MSSM 



3 

n 

i=l L 



OiiMs 
ai{M, 



R 



and = n 



1 L«.(M 



G 



(38) 



Here 6j = 3^ for i = 1y,2l,3c are the well known MSSM /3-function coefficients, 

bj = (13,2,4,-2) for j = l^.j^, 2^, 2^, 3c are the /3-function coefficients for the SUSYLR 
model |5|, and jiS are the anomalous dimensions for the LLLL operator, calculated in 
Appendix B. Using these results, we obtain 



A 



MSSM 



OsiMs 
asiM, 



R 



-4/3 



R 



aiAMs) 
aiy{M, 



R 



1/33 



0.91 



(39) 



using the MSSM running of the gauge couplings, and similarly. 



A 



SUSYLR 



'asiMR) 


-2 


«2JMk)' 


3/2 


«2H(Mfi)' 


3/4 








«2,(Mg)_ 




[«2«(Mg)J 




[a,,_,{MG)\ 



1/26 



0.08 (40) 



using the SUSYLR running of the gauge couplings jsl. Combining all these results, we get 
the overall extrapolation factor in bringing the operators from the GUT scale down to 1 
GeV: 



A = y^NS^MSSM^SUSYLR ^ Q H 



(41) 



D. Predictions for partial lifetimes 

Substituting the extrapolation factor obtained in Eq. fHTj) in the expression for the partial 
decay width given by Eq. and using Mt — Mjj ~ 4 x 10^^ GeV in our model, we obtain 
the partial lifetimes of different decay modes: 

_ ^ ^ (4.42 X 1033 years) / lO^^X / 200 GeV\ ' ( ^\ ' 
^"r,- |/(F,D)P [\C\^)[ j llTevj ^^'^ 
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The wino mass, m^y, has been constrained at LEP to be larger than ~ 100 GeV 31| . 
essentially independent of any specific model. As a typical value, we choose the universal 
gaugino mass, mi/2 = 200 GeV, which when extrapolated to the weak scale gives ~ 134 
GeV for the wino mass. 



Model (A) 

As we are interested in obtaining an upper bound on the partial lifetimes of various proton 
decay modes, we adopt the strategy of varying the mixing parameters Xj's defined by Eq. ( 128|) 
to maxirnize the expression f H2]) and simultaneously satisfying the present experimental lower 
bounds [32]. We find that the most stringent constraint comes from the p — > K^V decay 
mode, and for this decay rate to be consistent with the present experimental bound, we 
must have the sfermion mass > 1.2 (2.1) TeV for the MSSM tan /3 = 10 (30). This value 
of Mp when extrapolated to the GUT-scale, puts a lower limit on the universal squark mass 
mo for a given value of mi/2- The allowed region in the mo — mi/2 plane satisfying the proton 
decay constraints and also satisfying the EWSB constraints is shown in Fig. [5l It is clear 
that this model favors low values of tan/3. 

2400 
2200 
2000 
g 1800 

o 

^ 1600 
1400 
1200 
1000 

200 300 400 500 600 700 800 900 1000 
mi/2 [GeV] 




tan(3=10 ♦ 
tanp=30 



FIG. 5. Model (A) allowed region in the mo — rni/2 plane satisfying the proton decay and EWSB 
constraints for tan/3 = 10 (red) and tan/3 = 30 (green). 

The model predictions for the upper bound on partial lifetime of various proton decay 
modes are given in Table [Tll We also list the present experimental lower bounds for compar- 
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Decay 
mode 


Experimental 
lower limit (xlO^'^ yr) 


Predicted upper limit ( x 10'^^ yr) 


tan /3 = 10 


tan /? = 30 


p — )- K^V 


2.3 


2.3 


2.3 


p K^n+ 


1.3 


399.3 


738.8 


p K°e+ 


1.0 


1.3 X 10^ 


49.7 


p — TT^e~^ 


10.1 


5.8 X 10^ 


230.0 


p — )• 7r°^+ 


6.6 


2.4 X 10^ 


1.3 X 10"^ 


p — )• vr"*"!/ 


0.025 


1.5 


0.8 



TABLE II. Model (A) predictions for the upper limits on the partial lifetimes of various proton 
decay modes in 50(10) with low scale SUSYLR for tan/3 = 10 and 30 for mi/2 = 200 GeV. We 
have chosen the value of the universal scalar mass niQ to be 1.2 (2.1) TeV for tan/3 = 10 (30) so 
that the p — )• K'^V constraint is just satisfied. The present experimental lower limits are also given 
for comparison. 



ison. As noted above, the most stringent constraint on the parameter space comes from the 
p — )■ K^u decay mode; this is due to the fact that the neutrino final states add incoherently 
for the three generations, and hence, the decay rate for the neutrino final states will be 
much larger compared to the rates of other decay modes due to the third generation Yukawa 
coupling dominance. This also explains why the p — )■ tt^Tj decay rate is so large, even though 
it is Cabibbo-suppressed. The predicted upper bounds for these neutrino final states may be 
testable in the future proton decay searches, as in the next round of Super-Kamiokande [32 1 
or megaton type detector searches. 



Model (B) 

As in the model (A), we maximize the function |C|~^ given by Eq. (129]) with respect to 
the Xi parameters to find an upper bound on the proton decay lifetime. However, due to 
the particular structure of the Yukawa matrices in this model, as given by Eqs. (HM and 
( !T5|) . the parameters X2 and X3 have no effect on the amplitude and the only effective mixing 
parameter is Xi. The experimental lower bounds on the lifetime of various proton decay 
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modes will then put a lower bound on the ratio — —. It turns out that the most stringent 
bound is p — )■ K^v {n^^'^) for tan /3 = 10 (30) and we must have 

^ > 1.44 (1.06) X 10^ GeV (43) 

As an example, for mi/2 = 200 GeV and Xi = 0.1, it puts a lower bound on the first and 
second generation squark masses to be M? > 1.4 (1.2) TeV for tan/3 = 10 (30). The model 



Decay 
mode 


Experimental 
lower limit (xlO'^'^ yr) 


Predicted upper limit (xlO'^^ yr) 


tan 13 = 10 


tan /3 = 30 


p — )- K^V 


2.3 


2.3 


3.5 


p K^fi+ 


1.3 


2.3 


1.6 


p K^e+ 


1.0 


* 


* 


p — TT^e^ 


10.1 


* 


* 


p — )• 7r''/_i+ 


6.6 


9.8 


6.6 


P — )• TT^V 


0.025 


1.7 


2.7 



TABLE III. The predictions for the upper limits on the partial lifetimes of various proton decay 
modes for the new mass fit in our model for — 200 GeV and xi = 0.1. The most stringent 
constraint is from p — t- vr'^/U^ mode, and hence. The squark mass has been chosen to be 1.4 (1.2) 
TeV for tan /3 = 10 (30) so as to just satisfy the most stringent bound. Note that in this case, 
the model does not have any predictions for the decay modes p — t- K^e^ and p — t- vr'^e"'", because 
the C coefficients for both these modes involve products of (1,1) elements of the Yukawa coupling 
matrices, and by construction, these elements are zero for all the three coupling matrices; hence 
these modes have vanishing decay rates. 

predictions for xi = 0.1 for various decay modes are given in Table [ml We note that the 
observation of one of the decay modes in the last two columns of Table IIIII at a given rate 
will fix xi and the rates for remaining modes (the ones without stars) are then predicted and 
should provide a test of this model. It should also be noted here that within the mSUGRA 
framework at low tan /3, Tevatron has put a lower limit of 375 GeV for the squark mass 
based on an integrated luminosity of 1 fb~^. We expect our predicted lower bound on the 
squark mass which is of order 1 TeV to be testable at higher luminosities within the reach 
of LHC. 
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VII. EFFECT OF i?-PARITY BREAKING 



So far we assumed matter parity so that there is no i?-parity violating terms in the 
superpotential (i.e. W — 0). In this section we discuss the imphcations for relaxing this 
assumption on proton life time. This is an interesting exercise in view of the fact that in 
MSSM embedding into SU{5), relaxing i?-parity (or matter parity) conservation leads to 
new contributions to baryon number violation with arbitrary strength, so that in principle, 
such models are not viable without matter parity assumption. We would like to study in 
this section the situation in the case of our SO{10) model. 

The most general it!-parity violating interactions upto dimension-5 operators in our model 
are the following: 

Mpi 

where ipa,b,c denote matter spinors and ipH and ipH are Higgs spinor fields. Before proceeding 
to discuss their implications, note that must be of order TeV otherwise the right handed 
neutrino field would decouple from the low energy sector and break the gauge multiplet 
required to implement inverse seesaw. There are the following classes of i?-parity violating 
operators that follow from this in conjunction with the + Wsb at the TeV scale: 

W-'(TeV) = M'.Lir + AW + [QIQIQI + LaQbQl + ^L^Lc + • • •] (45) 

Mpi 

Note that the first three terms within the square bracket, after B — L breaking, give rise 
to the familiar MSSM i?-parity breaking terms with however couplings determined to be of 
order which is of order 10~^^. Hence their contribution to proton decay is negligible. 
Note this would not be the case with 5'O(10) models where B — L symmetry is broken at 
the GUT scale. 

VIII. CONCLUSION 

In summary, we have discussed proton decay as well as electroweak symmetry breaking in 
a new class of recently proposed SO{10) models with TeV scale Wr. We showed in an earlier 
paper that the model explains small neutrino masses via the inverse seesaw mechanism and 
has the feature of gauge coupling unification. The right-handed neutrinos in this model are 
almost Dirac type (pseudo-Dirac) with masses also in the TeV range making them (as well 
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as the Wr and Z' bosons) accessible at the Large Hadron Colhder. The coUider signals are 
different from the case with Majorana right handed neutrinos of conventional type I seesaw. 
We have explored two classes of fermion mass fits in these models. In both the cases, all the 
Yukawa couplings entering the dimension-5 proton decay operators are fixed within certain 
assumptions by charged fermion mass fits, thereby leading to definite expectations for the 
partial lifetimes of various proton decay modes. We find that it is possible to satisfy the 
current experimental lower limits on the lifetimes with a wino mass of 100-200 GeV and 
squark and slepton masses of order TeV. More specifically, to satisfy the most stringent 
bound coming from the p — )■ K'^V decay mode, we need to have a lower limit of 1.2 (2.1) 
TeV on the squark masses in the case of model (A) for tan/3 = 10(30) and similar lower 
bounds for model (B) for a given 10-Higgs mixing, assuming the universality of squark 
and slepton masses, as in a typical mSUGRA type scenario. Thus, discovery of squarks at 
LHC can throw light on the validity of these models. It is also worth pointing out that 
the choice of 5*0(10) multiplets in this class of models is derivable from fermionic string 
compactification. 
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Appendix A: RGEs for soft SUSY-breaking masses in SUSYLR model 

Assuming i?-parity conservation and the trilinear couplings A's and Y's in the superpo- 



tential and soft breaking Lagrangian given by Eqs. f[T6|) and f[T7j) to be zero, the soft breaking 
mass RGEs at one-loop level are given by 16[ 

1 32 1 

--MMqI - QM^lMIlqIl - ^M.Mlal + (Al) 



1 32 1 

-MiMlgl - QM2RMl^gl^ - —UMal - g^7i^2, (A2) 
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where 



,/3 



d 



dt 



d 



o 



IGtt — m 



dt 



167r^— m| 



(it 



a/3 



167r — m$ 



-3M,Ml9l - QM2RMl^gl^ + lglS2 

o 

4/i2:/x^c(mJ.+mi.)+2/i2V( 

,^ Tr + h'^hi) + Tr (s/^l/., + h'X) ^1., 

Tr (6/il/ibm|e + Qhlmlh + 2h'^himl + 2h'lmlh', 
{-QM2LMlLglL - 6M2RMlj,gl^) S^b 



ml 



+ 



S2 = 4 [Tr (ml - m^, -m\ + m^^) + (' 



m^c — m|, 



(A3) 



(A4) 

(A5) 
(A6) 
(A7) 



(A8) 



(A9) 



We have ignored the RG running of the couphng ^J^c as these are higher order effects. 



Appendix B: Anomalous dimensions of the dimension-5 operator 

Here we present the derivation of the anomalous dimensions of the dimension-5 operators 
of the LLLL type given by Eq. (125!) . The calculation is straightforward in a supersymmetric 
gauge due to the fact that the operator Ol is purely chiral (it is an F-term), and hence, 
it follows from non-renormalization theorems that in a supersymmetric gauge, it will only 
have wave function renormalization. Then it is easy to show that the anomalous dimensions 
of any purely chiral operator are given by 



lO = I]C'2(? 



(Bl) 



where C2{r) is the eigenvalue of the quadratic Casimir operator in the representation r, and 
the sum runs over all the chiral superfields occurring in the chiral coupling. As the gauge 
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bosons belong to the adjoint representation, we have 

C2{r 



2N 



for SU{N) 

X 



for ^7(1 



(B2) 



Thus we have for SU(?>)c 



73, = 3 X - = 4 



(B3) 



as there are three 5'f/(3)c fields in the LLLL operator [e.g. {qq){qiy\- Similarly, we have 



72i,« 
7ly 

7lfl-L 



4x- = 3, 

4 
1 
4 



3 _ 1 

5 ~ 5' 

3 _ 1 

2 ~ 2 



(B4) 
(B5) 

(B6) 



Here the factors | and | are the GUT normalization factors for U{1)y and U{1)b~l respec- 
tively. 

We note that the same results would have been obtained in a non-supersymmetric gauge, 
though the calculation is much more involved. For instance, the same results were obtained 
to.t.eMSSMc.e.n.We.-Zu.,„ogau.e,„Ref.Q. 



Appendix C: The hadronic factors f{F,D) 



As noted in Sec. VI, the hadronic factor f{F, D) estimates the chiral symmetry breaking 
effects on different final states. The low-energy parameters D and F are usually chosen to 



33j. Then D + F 



be the same as the analogous parameters in weak semileptonic decays 
g(rip) _ ]^_27 is the nucleon axial charge, while D — F = g^^ = 0.33 — 0.34 [29]. This gives 
D = 0.8 and F = 0.47. Using these constants and the approximations rriu^d ^ f^s ^ as 
well as — ^ where is the momentum transfer (the momentum of the anti-lepton 
for physical decays), all the hadronic matrix elements can be obtained 28|. In Table IV, we 
list the results for different decay channels. 
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